Vortex macroscopic superpositions in ultracold bosons in a double-well potential 
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We study macroscopic superpositions in the orbital rather than the spatial degrees of freedom, 
in a three-dimensional double-well system. We show that the ensuing dynamics of N interacting 
excited ultracold bosons, which in general requires at least eight single-particle modes and 0(N 7 ) 
Fock states for large N, is described by a surprisingly small set of many-body states. An initial state 
with half the atoms in each well, and purposely excited in one of them, gives rise to the tunneling of 
axisymmetric and transverse vortex structures. This process generates macroscopic superpositions 
only distinguishable by their orbital properties and within experimentally realistic times. 
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Orbital physics play a crucial role in many important 
phenomena, due to the combination of orbital degeneracy 
and anisotropy of the vibrational states, and its correla- 
tion to other attributes, like charge or spin. For example, 
in metal transition oxides it is key to understanding high- 
temperature superconductivity or colossal magnetorre- 
sistance [2[. Recent experiments paved the way for the 
exploration of their role in the physics of ultracold atoms 
in three dimensional (3D) optical lattices [3|, where this 
degree of freedom can be separated from those of charge 
and spin, and is the origin of exciting properties, such as 
novel phases or supersolidity [4]. Despite the fact that 
ultracold atoms are a natural system for realizing macro- 
scopic superposition (MS) states [l(, these have not been 
experimentally demonstrated therein, in part due to the 
very short decoherence times for such states. In this Let- 
ter we propose a new kind of MS state with the potential 
for greatly increased decoherence times, a vortex macro- 
scopic superposition (VMS), based purely on orbital an- 
gular momentum instead of spatial superpositions in a 
3D double well (DW). 

Ultracold bosons in DWs are conventionally described 
by a two mode approach, accounting for the atoms con- 
densed in the ground state and localized in either one of 
the two wells. If a great part of the population is inten- 
tionally excited to the first level of energies [3| , the three 
degenerate orbital states localized at each well and show- 
ing z-component of the angular moment m = —1,0,1 
have to be considered, together with new relevant pro- 
cesses, as sketched in Fig. [1] and discussed below. We 
assume that initially the atoms are distributed equally 
between both wells, and excited to an orbital state in one 
of them. Even though this eight-mode picture requires 
a great number of Fock states, we show that its dynam- 
ics is described by only a few many-body states. Then, 
while the number occupation of both wells remains con- 
stant under time evolution, the vortex tunnels between 
wells with a period shorter than the lifetime of a conven- 
tional experiment. This process is accompanied by the 
creation of VMSs. Spatial MSs decohere after a single 
interaction with an external particle. In contrast, VMSs 
must interact with many particles in order to spatially re- 
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FIG. 1. (Color online) Schematic of the 3D DW. The eight or- 
bital modes are represented as distorted spherical harmonics 
(magenta) and the processes among them by arrows. In the 
right well, the arrows represent both tunneling and interact- 
ing processes between different orbit als. The blue "peanut" 
shows the actual shape of an isoenergy surface of the DW. 

solve the presence of the vortex on one side or the other, 
as external particles such as are found in the imperfect 
vacuum of experiments, as well as electromagnetic fields 
used to create the DW do not have a well-defined angular 
momentum; thus VMSs are stronger against decoherence 
than other MSs in DWs Q. 

Ultracold bosons condensed in the ground state in 
DWs undergo two major processes: they interact in pairs 
in the same well with energy U or they can tunnel to 
the other well with energy J. The phenomena predicted 
by semiclassical mean field approaches .e.g., macroscopic 
quantum tunneling and self-trapping [6[, were observed 
already in experiments [7]. Also, it has been shown that 
vortices can tunnel in DWs and that vortex-antivortex 
MSs can be obtained in a single trap [8]. Moreover, 
topological defects have attracted recently attention in 
connection to MSs [9j, which in DWs are expected in 
the high barrier or Fock regime NU ^> U ^> J. In the 
Fock regime, mean field approaches cease to be useful, 
and other methods are required, with multiconfigura- 
tional Hartree methods the most numerically accurate 
techniques known to date [10], but impractical in the 3D 



DW for VMSs. A simpler method based on direct diago- 
nalization of Lipkin-Meshkov-Glick (LMG) like Hamilto- 
nians with more than two modes is used in this Article, 
because it permits one to treat the 3D DW, to calculate 
with more atoms, and to obtain analytical results with 
perturbation techniques 



U|, as we will show. 



Hamiltonian - A dilute gas of N ultracold bosons 
of mass M interacting through two-body interactions 
with coupling constant g = 47rfi 2 a s /M, where a s is 
the s-wave scattering length, is trapped in a 3D DW, 
which we assume is given by a separable potential V(x), 
x = (x,y,z) G M 3 , with harmonic ID potentials in the 
x and y directions, and a Duffing potential in the z di- 
rection, V{z) = V (Sz 2 /a 2 + 16z 4 /a 4 + 1). Since V(x) 
takes the shape of a peanut with lobes in the z direction, 
as shown in Fig. [TJ we call it the peanut potential We 
consider a two-level eight single-particle mode expansion 



of the field operator in second quantization [12j , giving 



H = X^m Him + H[ nt , with £ = 0, 1 the level index and 
AE = Ei—Eq the energy difference between levels. Here, 
Him are LMG Hamiltonians for each level: 
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where h] lm 
well index j 



°j£m L 



is the number operator with 



1,2 for the left and right wells, respec- 
tively, level £, and z component of the angular momen- 
tum ra. Here b*- irn and bj£ m are the creation/annihilation 
operators satisfying the usual bosonic commutation re- 
lations [l2j. The interaction and tunneling coefficients 
for atoms with the same £ and m are Ug m and J^ m , 
respectively (see Fig. [TJ. On the other hand H in t = 



Tin 



i^intra is given by 
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Equation (|2j) with m = accounts for a zero-vorticity 
interlevel hopping process that permits two atoms in the 
ground state to be excited to an orbital with m = 0, or 
vice versa, with an energy Uqq. For m = ±1, Eq. (J2j) ac- 
counts for a vortex- antivort ex interlevel annihilation and 
creation process, in which two atoms are excited to (or 
decay from) an orbital with m ^ 0, each with a different 
sign of m, with energy U 11 . Equation (J3]) accounts for 
a third process, vortex- antivort ex intralevel annihilation 
and creation, in which two excited atoms with m = can 
generate a pair of atoms with m = ±1 (each with differ- 
ent sign of m), or vice versa, with energy U 11 (see Fig. [I]). 
Notice that the 3D DW shows cylindrical symmetry with 
respect to the z axis, and two- fold (Z 2 ) symmetry with 



respect to the x and y axis. Thus the z component of 
the angular momentum, h z = J2mnji m , is conserved. 
Conversely, the angular momentum operator L 2 , whose 
lengthy expression can be obtained for 3D DW using a 
conventional procedure [13|, is not generally conserved. 

We focus on the Fock regime, for which the inter- 
action coefficients are much bigger than the tunnel- 
ing ones, Q m = Je m /Ue m <^ 1? and we further as- 
sume that x — ^ m ax^ m £/ m /[/7/ r ^ n }/AE <C 1, as is 
necessary for the eight mode single-particle approxima- 
tion to be accurate [12j. We expand our states ip in 



terms of the Fock basis, \ip) = ^Z i=0 Ci\i), where Q = 
[(N + 7)!]/[(iV!)(7!)] ~ N 7 /7\ for N > 1 is the di- 



mension of the Fock space and 
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\n-lJ = Kl)- 1/2 (^ m ) U ^|0). We focus on the dy- 
namics of vortex structures for zero population imbal- 
ance, i.e., the same number of atoms in each well. If all 
atoms remain in the lowest orbital, no tunneling process 
occurs. Conversely, if the atoms in one of the wells are 
excited to an orbital with the experimental techniques 
available [14|, vortex structures tunnel between wells, ac- 
companied with the creation of VMS states, as discussed 
below. There are two cases: (i) the vortex is axisymmet- 
ric and lies in the z-direction, i.e., m = ±1; and (ii) the 
vortex is transverse, lies in the x — y plane, and m = 0. 
Other limiting cases, including non-zero population im- 
balance, are discussed in [13j . 

Case (i): Axisymmetric Vortex - The initial state is 

\i) with n[% = N/2 and 4oo = N / 2 - We use hi £ h or " 
der perturbation theory to analyze the spectra of the 
two-level Hamiltonian, with the perturbing Hamiltonian 
the hopping terms in Eq. (pQ) , and the interlevel Hamilto- 
nian, Eqs. (J2j) and (J3j). We obtain quasidegenerate paired 
eigenvectors, ip = (1/V2)(\i) ± |j)), with n^ = N/2 and 
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TV/ 2, and splitting given by 



Ae ±1 = 2(N/2\fC 2 J^U, 
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U(n,p) = (-l)(»+*> {/7 ° °2n[7V/2 - n) + U^2p[N/2-p] 

-4U^[N/2(n+p)-2np}}-\ 

with N/2-1 < E^i 1 ij < N/2. We validated Eq. (g]) 
for small N in simulations, and note that, despite the 
complex sums and lack of a simple expansion, the es- 
sential scaling is Ae±i - U§§(JooJii) N/2 /Nl(U§§) N , 
since Uqq ~ U^ ~ U 11 up to constant factors of or- 
der unity. In this case, the atoms slosh between wells 
with period T = 27rfi/Ae±i, half of them always ex- 
cited with m = 1, the other half non-excited, with oc- 
cupation nj£ m (t) = ('0(t)|n J ^ m | / 0(t)) periodic for fim 
and n200- The numerical calculations for these two vari- 
ables via exact diagonalization are shown in Fig. [2fa), 




FIG. 2. (Color online) Axisymmetric Vortex Tunneling, (a) 
Average occupation of left well for unexcited, zero orbital an- 
gular momentum atoms (solid curve), and excited atoms in 
an axisymmetric vortex with m — 1 (dash-dotted curve), (b) 
Probability densities Pi (£) showing the dominant Fock states 
i contributing to time evolution. The excited and unexcited 
atoms slosh between wells with the same period, (c) Spatial 
(solid thin curve) and angular momentum (solid thick curve) 
entanglement entropy, (d) Total L 2 (solid curve) and total L 2 
(dash-dotted curve) are conserved, while on-site h z (dashed 
curve) is not. 

for Z7 ° = 1, Coo = 1/5, Jio/Joo = 5, x = 1/10. The 
probability P\ii(t) of finding N/2 excited atoms in well 
j = 1 with m = 1 and the probability P20o(^) of them 
being non-excited in the other well are equal and given 
by Pni(t) = P200O) = cos 2 (Ae 1 t/2h). At the quarter 
period Pm(t) = P2oo(t) = 1/2, and we obtain a VMS in 
the z-component of the angular momentum. Figure Efb) 
plots the probability density in time, |q(£)| 2 , showing 
that at t = T/4 the system is in a VMS. Since we consider 
only pure states, the VMS has zero total quantum en- 
tropy. However, the local entanglement or von Neumann 
entropy in space or angular momentum are both non- 
zero: the partial trace over the four modes in the right 
well, or alternately, over all modes but £ = l,m = ±1, 
yields algebraically complex expressions for Sj=i and 
S m =±i, which we don't show here for brevity [13]. We 
normalize to the maximal entropy in each case and plot 
the results in Fig. [2fc). Both show a maximum at the 
quarter period, corresponding to the dynamical occur- 
rence of the VMS. Finally, both total L 2 and total ~h z 
are conserved in this case. Conversely, on-site h z , ob- 
tained after restricting the index j to one of the wells, is 
not conserved and oscillates with period T. 

Case (ii): Transverse Vortex - If the initial state \i) 

has 77,110 = N/2 an d r4oo = N/2, with TV assumed even, 
the intralevel vortex-antivortex annihilation and creation 
process plays a crucial role, and pairs of atoms with m = 
±1 can be created from pairs of excited atoms with m = 
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FIG. 3. (Color online) Transverse vortex tunneling. Panel 
layout as in Fig. [2] (a) Average occupation of left well for un- 
excited atoms (solid curve), and atoms in transverse vortex 
(dashed curve), (b) Four many-body modes now dominate, 
rather than just two - compare to Fig. [2] (c) Spatial (solid 
thin curve) and angular momentum (solid thick) entangle- 
ment entropy, here for m — 0. (d) Both total L z (dashed 
curve) and on-site L z (superposed) are conserved, while total 
L 2 (solid curve) is not. 

0. Using perturbation theory we find that the relevant 
Fock vectors include not only \i) and |j), but also, due to 
this process, two more vectors \k) and \l), where n 



(k) _ 



N/2 and n% = N/2; < 00 = N/2, n™ = N/2 - 2, 



M _ 



100 



and 77-21 ±1 = 1' an d with tt^qq 



N/2, n% = N/2 - 2, 



and n[{ ±1 = 1. Thus despite 0(N 7 ) Fock vectors, the 
dynamics is described by combinations of just four of 
them, the quasidegenerate pairs ip± = a\i) ±a\j) + f3\k) ± 
P\l), with splitting Ae , and <j>± = =F/3|i) - j3\j) ± a\k) + 
a\l), with splitting Ae f . Perturbation theory shows that 
a ^> /3 and all other couplings to other Fock vectors are 
negligible [13]. We find the average number of particles 
in the right well in a transverse vortex state to be 

(n 20 o) = 2C(1 - cos Ae t) + 2C"(1 - cos Ae' t) , (5) 

where h = 1, C = a 2 N(a 2 + f3 2 ), C = f3 2 N(a 2 + f3 2 ), 
and C ^> C since a ^> j3. In Fig. [3] we present, for this 
case, (a) the evolution of the average occupation of well 
j = 1 for excited atoms either with m = or m = 1, 
and (b) of the probability density Pi(t), when /Too = 1, 
Coo = 10" 2 , J10/J00 = 5, TV = 8, and x = 1/10. This 
four eigenvector problem leads to a quasiperiodic motion 
where the two relevant frequencies are proportional to the 
splittings, and where the modulation is very small due to 
the small coupling to Fock vectors with riji±i = 1. The 
local entanglement entropies Sj=i and S m =o, obtained 
from the partial trace over all modes but £ = l,m = 0, 
are shown in Fig. [3] (c). Now, the atoms tunnel between 
both wells with the same period, creating a transverse 
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FIG. 4. (Color online) Energy splittings for VMSs. Axisym- 
metric (analytical: solid blue curve; numerical: blue crosses) 
and transverse (analytical: dashed, dotted, and dash-dotted 
red curves; numerical: red asterisks and circles) vortex cases. 
In the transverse case there are two limits: small and large 
1/C- We represent the numerical results for x — 1/80 (aster- 
iks) and x — 1/10 (circles). The analytical small l/£ limit 
coincides in both cases (dashed curve). The analytical large 
1/C limit depends on AE and is represented as a dotted (dash- 
dotted) curve for x = 1/80 (x = 1/10). 

VMS, and vortex- ant ivortex pairs are rapidly created and 
annihilated during this process. Finally, Fig. [3] (d) shows 
that, while L 2 is not constant, both total and on-site h z 
are conserved, since the atoms with m = ±1 are created 
in pairs. 

The splittings can again be obtained using perturba- 
tion theory. There are two possible scenarios. First, if 
Coo is bigger than x, the splitting will be given by Eq. ^ 
upon substitution of Jn by Jio, U\\ by U\q, and Uqq by 
Uqq. Second, if Coo is much smaller than x, the coupling 
between \i) and \j) is due to the zero-vorticity inter level 
hopping. Then, this splitting is 



io\^/2f r 



(6) 



A4 = 2(7V/2)(7V/2!) 2 (/7 1 °) iV/z /7, 

U(n,p) = {U^[4n(N/2 - n) - 2n - 2p(2p - 1)] + 

Ul°[4p(N/2-p-±)-2n(2n-l)]+2(p-n)AE}-\ 

with U given in Eq. ((4]), ij running only to j = 
N/2 — 1. Now, the essential scaling is A^o ~ 
(N/2)\(U^) N ^ 2 /N(N - l)AE N / 4 (U^) N ^-\ In Fig. H 
we show the splitting, both for the axisymmetric and 
transverse vortices, for TV = 8 and different values of Coo- 
The axisymmetric splitting is smaller than the transverse 
splitting, and the analytical approach shows good agree- 
ment with the numerical calculation: thus transverse vor- 
tices tunnel faster; moreover, a transverse VMS is easier 
to make in experiments [14|. We also plot the analyti- 
cal calculation for the splitting for the transverse vortex. 
First we use Eq. ((4]) as discussed above. The numer- 
ical calculation shows good agreement with this curve 
for small 1/Coo- For bigger values of this coefficient, the 
splitting tends to an asymptotic value, given by Eq. ((6]), 
which is increased for smaller values of \- 

Let us obtain the period of oscillation for typical ex- 



periments with uj = 2tt x 70 Hz to 7 Khz with AE = 
ujH. Taking x = V 2 we obtain *7 ° = x^E/N = 
(0.125tol2.5)ft Khz, which for Coo = 1/100 gives J 00 = 
(1.25tol25)ft Hz. For these quantities we obtain Ae w 
7 (10 -3 to 10 _1 )fi kHz, which gives a period of oscillation 
T = 2irh/ 'Ae = 1 to 0.01 s, that is the oscillations will 
occur with a frequency of 1 to 100 Hz. MSs will be ob- 
servable in an experiment if this time is shorter than the 
decoherence time. Conventional NOON states, where all 
atoms occupy simultaneously both wells, are extremely 
fragile against decoherence processes, such as those in- 
duced by imperfections of the potential [ll| , spontaneous 
emission, or the thermal cloud, among others [5[, since 
they decohere after a single interaction. Since these 
sources of decoherence do not have a well-defined orbital 
angular momentum, VMS states have a much longer de- 
coherence time: the core of the condensate has to be 
resolved by the environment to make the VMS collapse. 
The number of interactions is proportional to the total 
size of the condensate over the size of the core. For ex- 
ample, when the core occupies 1% of the total volume of 
the condensate, VMS decoherence time will be 100 times 
larger than for conventional MSs. 

In conclusion, we have shown how vortex macroscopic 
superposition states arise from time evolution of vortex 
tunneling in a double well potential. Initial homogeneous 
distributions of atoms in both wells, excited in one of 
them to an orbital with m = ±1 (axisymmetric vor- 
tex) or with m = ( transverse vortex), give rise to the 
tunneling of the vortex structure, while the total aver- 
age density remains constant in both wells. These VMS 
states are only distinguishable by their angular proper- 
ties, and thus represent a new route for the realization of 
MSs with ultracold atoms in DW with potentially much 
longer decoherence times. 
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